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Goals 


The  long-term  goals  of  our  coastal  research  at  Rutgers  include:  (1)  the 
implementation  and  testing  of  numerical  models  suitable  for  regional  problems 
featuring  irregular  coastlines,  steep  topography,  and  strong  coastal/deep  ocean 
interaction;  (2)  the  use  of  these  models  to  simulate  the  regional  response  of  coastal 
currents,  fronts,  and  eddies  in  the  California  Current  System;  (3)  the  simulation  of 
storm-driven  sediment  transport  and  wind-driven  upwelling  in  the  Middle  Atlantic 
Bight;  (  4)  the  determination  of  the  local  balances  of  heat,  mass  and  energy  in  these 
regions,  including  the  patterns  and  methods  of  their  exchange  across  the  shelf/deep 
ocean  boundary;  and  (5)  the  application  of  these  regional  models  to  the  study 
of  coupled  physical/biogeochemical  processes.  ONR  N00014-91-J-1557  produced 
results  pertinent  to  issues  (1),  (2),  (4)  and  (5). 


Results  Under  ONR  N00014-91-J-1557 

Using  the  semi-spectral  primitive  equation  model  of  Haidvogel  et  al.  (1991b; 
development  supported  by  N00014-86-K-0751),  the  evolution  of  a  forced,  surface- 
intensified  eastern  boundary  current  has  been  studied  in  the  presence  of  both  finite- 
amplitude  topography  and  irregular  coastal  geometry  (Haidvogel  et  ah,  1991a). 
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The  laminar  surface  current  quickly  evolves  away  from  its  smooth  intial  conditions 
and  approaches  a  turbulent,  time-dependent  eciuilibrium  featuring  an  intense, 
meandering  along-shore  jet  with  local  velocities  of  80-100  cm/sec.  A  deep,  poleward 
undercurrent  also  forms.  Subsequent  interaction  with  the  protruding  cape  geometry 
causes  an  offshore  deflection  in  the  steepening  frontal  meanders,  some  of  which 
produce  elongated  filaments  which  penetrate  significant  distances  (400-500  km) 
offshore.  The  emerging  filaments  are  characterized  by  a  strong  downwelling  signal 
(maximum  vertical  velocities  of  30-40  meters  per  day).  The  simulated  filaments 
ultimately  pinch  off.  typically  within  40-50  days,  to  form  a  co-rotating  pair  of 
detached  eddies.  The  existence  of  the  cape  geometry,  as  well  as  the  southward 
surface  flow,  are  necessary  for  filament  generation:  removal  of  the  irregular  coastline 
or  reversal  of  the  sense  of  the  surface  circulation  inhibit  filament  formation.  Detailed 
analysis  of  the  instantaneous  and  time-mean  momentum  balances  shows  strong 
dynamical  similarity  to  observations  taken  during  the  Coastal  Transition  Zone 
experiment,  and  confirms  that  eddy  transport  mechanisms  are  responsible  for  the 
maintenance  of  the  poleward  undercurrent.  Offshore  transport  of  heat  by  the 
filaments,  found  to  be  0(0.02)  petawatts,  is  substantial. 


General  transport  patterns  and  residence  times  in  the  upper  100  meters  of  the 
Coastal  Transition  Zone  (CTZ)  were  studied  using  simulated  Lagrangian  drifter 
experiments  (Hofmann  et  ah.  1991).  The  circulation  fields  used  for  the  experiments 
were  obtained  from  numerical  simulation  described  above.  The  composite  drifter 
trajectories  show  that  transport  patterns  and  residence  times  in  the  upper  100 
meters  are  determined  by  the  proximity  of  the  drifter  release  point  to  the  offshore¬ 
flowing  filament.  The  magnitude  of  the  offshore  displacement  of  the  simulated 
drifters  implies  substantial  cross-jet  exchange.  In  particular,  comparison  of  the 
simulated  drifter  trajectories  with  the  trajectory  followed  by  a  drifter  released 


during  the  1988  CTZ  field  experiment  supports  the  hypothesis  that  the  deepening  ™ 
of  the  abundance  maxima  observed  in  certain  zooplankton  distributions  resulted 
from  these  populations  being  downwelled  as  they  were  transported  offshore  by  the 
flow  associated  with  the  filament  observed  in  1988.  ^  .  /)  ' 
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The  three-dimensional  CTZ  model  has  been  coupled  to  a  nine-component  food 
web  model  and  a  bio-optical  model  (Moisan  et  ah.  1995).  The  resulting  simulated 
rhree-ciimensional  nutrient,  plankton  and  submarine  light  fields  agree  well  with 
those  observed  within  the  CTZ.  Specifically,  high  nutrient  and  plankton  biomass 
occur  onshore  and  within  the  core  of  the  simulated  filament.  The  depth  of  the  1 
percent  light  level  shallows  to  less  than  30  meters  in  regions  of  high  phytoplankton 
Itiomass  and  deepens  to  greater  than  75  meters  in  regions  of  low  phytoplankton 
biomass.  The  integrated  carbon  transport  for  the  model  varied  with  distance  from 
shore  and  is  highest  in  the  region  where  the  filament  circulation  pattern  develops  a 
counter-rotating  pair  of  offshore  eddies.  The  annual  integrated  carbon  transport  is 
estimated  to  be  1.89  x  10^^  grams  C. 

The  model  used  by  Haidvogel  et  al  (1991a)  to  simulate  filament  evolution  in 
the  Coastal  Transition  Zone  has  been  redesigned  to  permit  accurate  simulation  of 
coupled  coastal/deep  ocean  interaction  on  both  regional  and  basin  scales.  A  semi- 
implicit  coastal  ocean  circulation  model  using  a  generalized  topography-following 
coordinate  system  has  been  implemented,  shown  to  perform  well  on  representative 
test  problems,  including  the  evolution  of  the  diurnal  mixed  layer.  These  tests, 
conducted  in  both  one  and  two  dimensions,  show  that  the  model  handles  mixed 
layer  cycling  quite  well  even  in  the  presence  of  strong  variations  in  topography 
(maximum  to  minimum  depth  ratios  of  order  100;  Song  and  Haidvogel,  1994a).  A 
successive  set  of  simulations  of  the  California  Current  System  featuring  alongshore 
variations  in  wind  stress  and  coastline  geometry  have  also  been  conducted  (Song 
and  Haidvogel,  1994b). 


Publications 

Haidvogel,  D.B.,  A.  Beckmann  and  K.S.  Hedstrom,  1991a;  Dynamical  simulations 
of  filament  formation  and  evolution  in  the  Coastal  Transition  Zone.  J.  Geophys. 
Res.,  96,  15,017-15,040. 

Haidvogel,  D.B.,  J.  Wilkin  and  R.  Young,  1991b;  A  semi-spectral  ocean  circulation 
model  using  vertical  sigma  and  orthogonal  curvilinear  horizontal  coordinates. 
J.  Comp.  Phys.,  94,  151-185. 


3 


Hofmann.  E.E..  K.S.  Hedstrom.  J.R.  Moisan.  D.B.  Haidvogei  and  D.L.  Mackas. 
1991:  Use  of  simulated  drifter  tracks  to  investigate  general  transport  patterns 
and  residence  times  in  the  Coastal  Transition  Zone.  J.  Geophys.  Res..  96. 
15.041-15.052. 

Moisan.  J.R..  E.E.  Hofmann  and  D.B.  Haidvogel.  1995:  ivlodeling  nutrient  and 
plankton  processes  in  the  California  Coastal  Transition  Zone.  2.  A  three- 
'iiemnsional  physical-biological  model.  J.  Geophys.  Res.,  submitted. 

Song.  Y.  and  D.B.  Haidvogel,  1994a:  A  semi-implicit  ocean  circulation  model  using 
a  generalized  topography-following  coordinate  system.  J.  Comp.  Phys.,  115. 
228-244. 

Song.  Y.  and  D.B.  Haidvogel,  1994b:  Numerical  simulations  of  the  CCS  under  the 
joint  effects  of  coastal  geometry  and  surface  forcing.  In:  Estuarine  and  Costal 
Modeling  III.  Malcolm  Spaulding  et  ah.  eds.:  ASCE.  New  York;  682  pp. 


4 


Reprimed  from  Estuarine  and  Coastal  Modeling  III 
Proceedings  of  the  3rd  Imemattonai  Conference 
Sponsorea  by  the  Waterway,  Port.  Coastal  and  Ocean  Div./ASCE 
September  S-IO,  1993.  Oak  Brook.  Illinois 

Numerical  Simulations  of  the  CCS  under  the  Joint 
Effects  of  Coastal  Geometry  and  Surface  Forcing 


\’iilie  Song' 


Dale  HaidvogeF 


Abstrarf. 

II  S-Coordinate  Rutgers  University  Model  (SCRUM)  of  Song  and 

Haidvogei(1993),  a  coupled  coastal  (86  m  depth)  and  deep  ocean  (4600  m  depth) 
California  Coastal  region  has  been  investigated  for  the  dynamics  of  the  California 
Current  System  (CCS)  under  the  joint  effects  of  realistic  coastline  geometry  bot- 
tom  topography,  wind  forcing,  and  surface  heat  flux.  The  model  equations  are 
the  three-dimensional,  free  surface,  primitive  equations  with  orthogonal  curvi- 
inear  cwrdinates  in  the  horizontal  and  the  s-coordinate  in  the  vertical.  Three 
kinds  of  surfa^  forces-steady  wind  field  with  positive  curl,  long-term  oscilla¬ 
tory  wind,  and  diurnal  cycling  wind  and  thermal  forcing-are  used  to  simulate 
he  variations  of  the  CCS.  Two  alternate  parameterizations,  dynamic  instability 
mixing  and  turbulence  closure  mixing,  are  used  in  the  diurnal  cycling  problem 
The  model  is  shown  to  be  stable  and  capable  of  handling  the  joint  effects  of 
geometry  and  surface  forcing.  Our  simulations  show  that:  (1)  upweiiing  and 
a  depressed  surface  height  are  always  found  at  the  coast  for  equatorward  wind 
stress;  (2)  the  equatorward  surface  jet  and  the  poleward  undercurrent  persist 
or  dl  upwelling-favorable  winds;  (3)  a  meandering  jet  and  eddies  can  be  caused 
by  the  coastline  geometry  and  wind  curl  in  combination;  (4)  surface  and  bottom 
l<i>ers  mixing  increases  the  variations  of  the  CCS. 

Introduction 

There  is  a  great  challenge  n,  modeling  the  coastal  oceans  due  to  irregular 
coastline  geometrv  and  severe  bottom  topography,  the  latter  featuring  both 

o^  r  '"'T'  geometric  rom- 

|.lc.Mty.  he  coastal  oceans  arc  the  sites  of  a  variety  of  physical  and  coupled 

phtsical/hiologicai  phenomena.  For  example,  eastern  boundary  regions  are  of- 
en  sites  of  wmd  -driven  coastal  upweiiing  and  surface  height  slopes  downward 
lo^vald  the  coast,  and  strong  vertical  advection  and  turbulent  mixing  due  to  sur- 
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Use  of  Simulated  Drifter  Tracks  to  Investigate  General  Transport  Patterns 
and  Residence  Times  in  the  Coastal  Transition  Zone 

Eileen  E.  Hofmann.'  Katherine  S.  Hedstrom.--’  John  R.  Moisan.' 

Dale  B.  Haiovogel,' and  David  L.  .M.ackas' 

General  tran.sport  ptitterns  and  residence  limes  in  the  upper  100  m  of  the  coastal  transition  zone 
(CTZ)  were  studied  using  simulated  Lagrangian  drifter  expenments.  The  circulation  heids  used  for  the 
drifter  experiments  were  obtained  from  a  regional  primitive  equation  model  that  incorporates  coastal 
geometry  and  botiom  topography  that  is  representative  of  the  CTZ  region.  The  simuiaied  circulation 
(ieids  show  velocity  patterns  that  are  consistent  with  those  associated  with  the  hlamenis  and  jets 
observed  in  the  CTZ  region.  Drifters  were  released  at  «80  points  in  the  model  domain  in  a  pattern  that 
bracketed,  both  horizoniailv  and  vertically,  the  region  in  which  olfshore-tiowing  hiaments  were 
observed  to  form.  At  each  location,  drifters  were  released  at  depths  of  30,  60.  90  m  and  tracked  for  30 
days.  The  depth  of  these  dnftcrs  vaned  along  their  trajectories  in  response  to  the  venical  velocities 
experienced  by  the  drifters.  An  additional  set  of  drifters  was  released  at  30  m  and  was  constrained  to 
remain  at  this  depth.  The  composite  drifter  trajectories  show  that  transport  patterns  and  residence 
limes  in  the  upper  100  m  are  determined  by  the  proximity  of  the  drifter  release  point  to  the 
otfshore-flowing  filament.  Drifters  released  inshore  of  the  filament  are  transported  southward  with  the 
mean  flow  of  the  model  California  Current.  Drifters  released  near  the  region  in  which  the  filament 
originates  are  transported  rapidly  offshore.  The  point  at  which  this  difference  in  drifter  fate  occurs  is 
approximately  100  km  offshore.  The  time  required  for  a  drifter  to  be  transponed  from  the  nearshore 
region  to  the  outer  boundary  of  the  area  covered  in  the  CTZ  model  domain  is  approximately  20  days. 

Comparison  ol  the  free-floating  drifters  with  the  30-m  fixed^depth  drifters  shows  that  the  deeper 
drifters  are  displaced  southward  as  they  are  transported  offshore  in  the  jet.  The  magnitude  of  the 
southward  displacement  for  the  90-m  drifters  is  about  15  km,  implying  substantial  cross-jet  exchange. 

Comparison  of  the  simulated  drifter  trajectories  with  the  trajectory  followed  by  a  drifter  released 
during  the  1988  CTZ  field  sampling  program  support  the  hypothesis  that  the  deepening  of  the 
abundance  maxima  observed  in  the  distributions  of  Dolioletta  ^ef^enhanri.  juvenile  forms  of  Euphau^ 
sia  pacifica  and  Eucalamts  californicus  resulted  from  these  populations  being  downwelled  as  they 
were  transported  offshore  by  the  flow  associated  with  the  filament  observed  in  1988. 


L  Introduction 

The  coastal  transition  zone  (CTZ)  off  the  coast  of  Califor¬ 
nia  (Figure  1)  is  characterized  by  meanders  of  the  California 
Current  and  by  prominent,  cross-shelf  jets  and  filaments. 
The  cold  filaments  can  be  identified  in  Advanced  Very  High 
Resolution  Radiomelry  [Breaker  and  Gilliland.  1981:  Ikeda 
and  Emery.  1984)  and  Coastal  Zone  Color  Scanner  imagery 
[Abbott  and  Zion.  1985,  1987]  and  arc  observed  to  extend 
several  hundred  kilometers  offshore.  Topically,  the  cold 
filaments  arc  associated  with  offshorc-tlou  ing  jets  that  have 
maximum  speeds  of  approximately  50  cm  s"'.  Several 
hydrographic  surveys  [Rieneckcr  et  aL.  I9S5:  Fhunent  et  al.. 
1985;  Kosro  and  Mayer.  1986:  Rienecker  and  \fooers.  1 9891 
have  been  conducted  to  describe  the  structure  of  the  cold 
filaments  and  jets  observed  off  Point  .Arena.  These  studies 
show  that  the  vertical  extent  of  the  offshore  jet  reaches  to 
approximately  200  m,  the  offshore  jet  is  about  50  km  wide, 
and  often  onshore  flow  is  found  to  the  south  of  the  jet.  A 
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similar  structure  was  described  for  a  cold  filament  off  Cape 
Blanco,  Oregon  [Moam  et  al.,  1988].  The  various  conceptual 
models  that  have  been  proposed  for  the  formation  and 
maintenance  of  the  cross-shelf  jets  and  filaments  are  sum¬ 
marized  by  Strab  et  al.  [this  issue).  The  offshore-flowing  jets 
represent  a  potentially  important  mechanism  for  transport¬ 
ing  properties  such  as  nutrients  and  organic  carbon  from 
nearshore  to  offshore  regions.  Recently,  Bucklin  et  al.  [1989] 
and  Mackas  et  al.  [this  issue  |  have  suggested  that  the 
filaments  are  important  in  determining  the  transport  and 
structure  of  zooplankton  populations  off  northern  California. 

In  1987  and  1988.  a  multidisciplinary  study  of  the  cold 
filaments  and  jets  in  the  CTZ  was  undertaken  [Coastal 
Transition  Zone  Group.  1988].  These  studies  have  provided 
hydrographic  and  velocity  measurements,  as  well  as  obser¬ 
vations  of  the  distribution  of  nutrients,  phytoplankton,  and 
zooplankton  in  the  CTZ  study  region.  Many  of  the  results  of 
the  1987  and  1988  CTZ  program  are  summarized  by  Strut  et 
al.  [this  issue).  As  part  of  the  1988  CTZ  field  observations 
(June-July  1988),  a  series  of  drifters  were  released  from  late 
June  to  early  July  in  a  cross-shelf  jet  [cf.  Strut  et  al.,  this 
issue).  The  drifter  study  was  similar  to  that  done  in  the  CTZ 
region  in  1987  which  is  described  by  Paduan  and  Niiler 
[1990].  As  in  1987  [Attott  et  al.,  1990],  biological  and 
hydrographic  observations  were  made  at  frequent  intervals 
(i.e.,  1-2  days)  at  the  drifter  location  (Figure  2).  The  assump¬ 
tion  underlying  these  observations  was  that  the  drifter 
tracked  the  same  water  mass  during  the  period  encompassed 
by  the  sampling. 

In  addition  to  the  CTZ  observational  program,  there  has 
been  a  complementary  modeling  component.  As  part  of  the 
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Dynamical  Simulations  of  Filament  Formation  and  Evolution 
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Using  the  semispecirnl  pnmiiive  cquaiion  model  of  Haidvogel  el  al.  (1991).  the  evolution  of  a  forced, 
surfacc-intensihed.  eastern  boundary  current  is  studied  in  the  presence  of  both  finite-amplitude 
topography  and  irregular  coastline  geometry.  The  model  domain  is  KXX)  km  in  alongshore  length,  and 
extends  on  average  7tK)  km  in  the  cross-shelf  direction.  A  representative  cape,  as  well  as  smoothed 
continental  shelf-siope  topography,  arc  included.  The  model  is  forced  by  inclusion  of  nudging  terms  in 
(he  equations  of  motion  which  relax  the  fiuid  system  back  to  a  prescribed  reference  state  on  a  time 
scale  of  45  da\  s.  The  reference  state  chosen  is  a  broad,  gcostrophically  balanced,  equaiorward  flow 
having  a  maximum  current  at  the  surl'ace  of  0.45  ms*  and  a  transport  of  approximately  10  Sv.  No 
explicit  wind  forcing  is  included.  Initialized  with  the  smooth  surface  current,  the  model  quickly 
approaches  a  turbulent,  time-dependent  equilibrium  featuring  an  intense,  meandering  alongshore  jet 
with  local  velocities  of  0.8-1. 0  ms'*.  A  deep,  poleward  undercurrent  also  forms.  Subsequent 
interaction  with  the  protruding  cape  geometry  causes  an  offshore  deflection  in  the  steepening  frontal 
meanders,  some  of  which  produce  elongated  filaments  which  penetrate  significant  distances  (400-500 
km)  offshore.  The  emerging  filaments  are  characterized  by  a  strong  downweiling  signal  (maximum 
vertical  velocities  ol  30—40  m  d  "  * ).  The  simulated  filaments  ultimately  pinch  off,  typically  within  40-50 
days,  to  lorm  a  corotatmg  pair  of  detached  eddies.  The  existence  of  the  cape  geometry,  as  well  as  the 
southward  surface  flow,  appears  to  be  necessary  in  this  model  to  produce  filament  generation:  removal 
of  the  irregular  coastal  gcomv'try  or  reversal  of  the  sense  of  the  surface  circulation  is  shown  to  inhibit 
filament  tormaiion.  Detailed  analysis  of  both  instantaneous  and  time-mean  momentum  balances  show 
dynamical  similarity  to  observations  taken  during  the  Coastal  Transition  Zone  experiment  and 
elucidate  the  eddy  transport  mechanisms  responsible  for  the  formation  of  the  poleward  undercurrent 
in  these  experiments.  Offshore  transport  of  heat  by  the  filaments,  found  to  be  0(0.06)  PW,  is 
substantial. 


I.  Introduction 

Observational  evidence  from  the  Coastal  Transition  Zone 
(CTZ)  experiment,  together  with  that  from  prior  in  situ  and 
remotely  sensed  measurements  collected  during  the  Coastal 
Ocean  Dynamics  Experiment  (CODE).  Ocean  Prediction 
Through  Observations  Modeling  and  Analysis  (OPTOMA) 
program,  and  other  programs,  provide  for  the  first  lime  a 
detailed  description  of  the  cold,  nutricnl-rich,  olfshore- 
Howing  filaments  which  form  during  the  spring  and  summer 
in  the  California  Current  system  (CCS).  Most  vividly  seen  in 
high-resolution  satellite  images  of  sea  surface  temperature 
(e.g..  Figure  1),  the  filaments  form  near  the  coast,  often  in 
proximity  to  regions  of  strong  coastline  curvature,  and 
penetrate  hundreds  of  kilometers  across  the  shclf-slopc 
boundary  into  the  deep  ocean.  The  filaments,  although 
narrow  (0(50-75  km)  in  width)  and  trapped  within  a  few 
hundred  meters  of  the  surface,  are  nonetheless  as.sociaied 
with  intense  offshore  currents  which  may  approach  a  meter 
per  second.  Strong  evidence  for  systematic  downweiling  of 
water  parcels  by  the  filaments  also  exists.  The  offshore  and 
vertical  fluxes  of  mass,  heat,  and  biological  materials  asso¬ 
ciated  with  these  filaments  can  therefore  be  substantial  and 
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are  of  likely  importance  to  regional,  and  perhaps  global, 
physical  and  biological  balances. 

The  CTZ  program  was  carried  out  in  1987  and  1988  to 
further  refine  our  understanding  of  the  dynamical  character, 
temporal  variability,  and  physical-biological  implications  of 
the  filaments  [Strub  et  al..,  this  issue).  In  addition  to  satellite 
imagery,  extensive  hydrographic,  biological,  microstruc- 
lure,  and  moored  current  meter  measurements  were  col¬ 
lected.  (See,  for  example,  Kosro  et  al.  [this  issue]  for  a 
description  of  the  large-scale  survey  cruises  in  1987.)  The 
Lagrangian  structure  of  the  filaments  was  also  investigated 
using  surface  drifter  observations  [Brink  et  al..  this  issue). 
Two  examples  of  the  types  of  data  thus  sampled  are  sum¬ 
marized  in  Figures  2  and  3.  which  show  inferred  along-shore, 
gcoslrophic  currents  from  the  1988  survey  and  observed 
l.agrangian  flow  structure  during  July  1988,  respectively. 
Further  reference  to  the  observed  structure,  and  dynamical 
character,  of  the  filaments  is  also  made  below  in  section  6. 

Despite  the  intensive  observational  effort  represented  by 
the  CTZ  program,  the  essential  nature  of  the  filaments  is  still 
open  to  interpretation.  For  example,  Strub  et  al.  identify 
three  conceptual  models  which  might  account  for  aspects  of 
the  observed  filament  behavior.  These  are  the  squirt  model 
(a  one-way  jet,  transporting  water  of  coastal  origin  to  the 
deep  ocean),  the  eddy  model  (in  which  a  slowly  evolving  set 
of  mesoscale  eddies  extracts  fluid  from  the  coastal  region), 
and  the  meandering  jet  model  (in  which  steepening  meanders 
on  a  strong  southward  jet  extend  offshore  and  pinch  off, 
thereby  transporting  coastal  fluid  to  the  deep  ocean).  None 
of  these  conceptual  scenarios  is  ruled  out  by  the  present 
observational  data  base:  however,  both  the  1987  and  1988 
observational  surveys  are  in  principle  consistent  with  the 
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